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ABSTRACT 



Aims. We implemented a fortran code that determine fundamental parameters of solar type stars from a list of Fe line equivalent 
widths. The solution should verify 3 conditions in the standard method: ionization equilibrium, excitation equilibrium and indepen- 
dence between metallicity and equivalent widths. We added the condition that the input metallicity of the model atmosphere should 
be similar to the output metallicity derived with equivalent widths. 

Methods. Solar-scaled Kurucz model atmospheres with NEWODF opacities are calculated with an independent program. Parameter 
files control different details, such as the mixing-length parameter, the overshooting, the damping of the lines and the weight factors 
in the definition of the^ 2 function. 

Results. FUNDPAR derive the uncertainties following 2 methods: the criteria of Gonzalez & Vanture (1998) and the dispersion us- 
ing the x 1 function. The code use the 2009 version of the MOOG program. The results derived with FUNDPAR are in agreement 
with previous determinations in the literature. In particular we obtained the fundamental parameters of 58 exoplanet host stars. The 
program is freely available from the web* . 

Key words. Stars: fundamental parameters - Stars: abundances 



1. Introduction 

Different methods have been used in the literature to derive fun- 
damental parameters and metallicities of solar type stars. For in- 
stance, some studies begin with a photometric estimation of tem- 
perature and gravity and then derive the metallici ty using equiv- 
alent widths. Useful codes such as WIDTH9 (iKuruczl Il993l 
1 1995b or BLACKWEfJ] (an implementation o f the BlackweU 
metho d) are widely used in the literature (e.g. Saff e & Levatol 
120041) . However the temperature and gravity are ussually consid- 
ered as fixed parameters. The observed stellar spectra could be 
comp ared with a grid of previously calculated sy nthetic spectra 
(e.g. [Fischer & Valenti I l2005"t ISaffe et al. 1120081) to determine 
fundamental parameters. However in this case the instrumental 
broadening should be taken into account and the rotational ve- 
locity could b e considered a s anoth er independent parameter. 

Recently, Sou sa et al. ] d2oToh derived an effective tem- 
perature calibration based on line equivalent width ratios of 
different absorption lines. Also the equivalent widths of Fe 
lines could be used to determine the parameters of solar 
type stars. The solution should verify three conditions in 
the standard method: [FeI/H]=[FeII/H] (i.e. ionization equi- 
librium), independence of the metallicity with the excitation 
potential (i.e. excitation equilibrium) and with respect to the 
equivalent widths. This method have been applied to so- 
lar type stars, for instance, in the determination of metallic- 
ity of stars wi th and without low mass companions or ex- 



opla nets (e.g . IGonzalez 1 1 19971 
Santos et al.T EOOO: 



1998 



Gonzalez & Laws 



IGonzalez et al. I Il999t 
120001: IGonzalez et al. I 



2001: lLaws & Gonzalez1l200lt Laws et al. 1120031) . The process 
estimate initially the fundamental parameters (T e ff> l°gg, [Fe/H] 



and the micro turbul ence velocity ). This information is used in 
the ATLAS dKurucz I [l993l 1 19951) program to derive a model 
atmosphere in LTE (local thermodynamic equilibrium). The 
model atmosphere with the measured equivalent widths in the 
spectra, are introduced in the code MOOG0 to derive a new 
metallicity. If the mentioned conditions are not satisfied, the pro- 
cess is restarted using new fundamental parameters calculated 
with the downhill method. 

In this contribution, we present the fortran code FUNDPAR 
(and their complement atlas. launcher) that reproduce the method 
explained. The programs are available from the wetQ, including 
detailed installation instructions and some technical details such 
as the format of the input/output files (install.txt). As an example 
of practical use, we derived the fundamental parameters of 58 
main sequence exoplanet host stars and verified the metal-rich 
nature of the group. The values derived are in agreement with 
previous determinations from literature. 

In the section 2 we show the general idea and the logic of 
the program. The procedure of minimization of the x 2 function 
is detailed in the section 3. The estimation of the uncertainties in 
the parameters and the comparison with literature, are showed 
in the sections 4 and 5, respectively. Finally we present some 
concluding remarks in the section 6. 

2. The logic of the program 

The algorithm is organized in one main procedure and a num- 
ber of sub-programs with specific functions. The iterative pro- 
cess begins in one starting point in a 4D-space, where the vari- 
ables are (T e ff< logg, [Fe/H], f). This starting point could be 



http://www 1 .appstate.edu/dept/physics/spectrum/spectrum.html 



2 http://verdi.as.utexas. edu/moog.html 

3 http://icate-conicet.gob.ar/saffe/fundpar/ 
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optionally given in the first line of the input file (see the file 
install.txt for details of input/output files). If this point is un- 
known, the program adopt (5000 K, 4.00 dex, -0.10 dex, 1.00 
km/s) to begin the iteration. For example, Sant os et al. I (2000) 
estimate T e ff and logg f rom the uvby photometry and the cali- 
bration ofJ01seirL(Q98J)^lvne^ 

lowing ISchuster & Nissenl (1 19891) and lEdvarsson et al. I (1 19931) . 
respectively. 

Then for each iteration step, i.e. for each group of 4 parame- 
ters (T e ff, logg, [Fe/H], £), the program should: 

a) Generate an appropriate ATLAS model atmosphere (through 
the atlas. launcher program). First, the code should select from 
the Kurucz grid the closer model atmosphere in the parameter 
space to those requested and take this as initial input model. 
Then the program execute ATLAS9 to derive the final model. 

b) Transform the model atmosphere in a format readable 
by MOOG (subroutine kurucz2moog). The format of the model 
atmosphere used by MOOG is not exactly the Kurucz model 
and should be rewritten accordingly. 

c) Call the MOOG program. In this point, the program 
takes the file containing the equivalent widths of Fel and Fell 
lines of the star and the model atmosphere as input for the 
MOOG program. The MOOG program is executed using a 
driver called "abfind" which is selected for the abundance 
determination. 

d) Read the new metallicity values calculated by MOOG 
(subroutine rmr, read-moog-results) and finally, 

e) Determine the value of the x 2 function. This step will 
be explained below. 

We take into account the conditions mentioned in the intro- 
duction in a variable called^ 2 . We adopt for;^ 2 the expression 
X 2 = w\c 2 + wic\ + W3C 2 + W4C4, where Wi,...,W4 are consid- 
ered weight factors (w, >=0), Cj and C2 are the slopes in the 
plots of [Fe/H] vs. logio(W/A) (logarithm of the reduced equiv- 
alent width) and [Fe/H] vs. excitation potential, C3 = [FeI/H]- 
[Fell/H], and C4 is the difference between the input ATLAS 
metallicity (step a) and the resulting metallicity using equivalent 
widths (step d). We added explicity the fourth condition: the in- 
put metallicity of the model atmosphere should be similar to the 
output metallicity derived with equivalent widths i.e. the term 
with C4. Then, the 4 conditions are quantified in the^ 2 function: 
the solution correspond to the minimum value of^ 2 . 

The user is free to modify the values of the weights Wi,...,W4 
under their own criteria. However we show a brief example esti- 
mating aproximately the values of the weights. Adopting^ 2 =1 
as the limit case of a solution, each condition contribute, for ex- 
ample, with 0.25 to the sum^ 2 = 0.25+.. .+0.25. In this case the 
4 conditions are taken equally important within^ 2 , which is not 
always true. In the plot of abundance vs. excitation potential, we 
accept a maximum slope, for instance, of c\ ~ 0.015/4 dex/eV, 
taking a difference of ~0.015 dex in abundance for a total range 
of ~4 eV in the excitation potential of Fe lines. Then in the limit 
case, w\c\ ~ 0.25 and thus w\ ~ 18000 eV 2 /dex 2 . The units of 
wi are forced to obtain the product w\c\ without units. In the 
plot of abundance vs. logio(W//l) (where W and A are the equiv- 
alent width and wavelength in A, respectively), we accept a max- 
imum slope (for example) of C2 ~ 0.015/1.5, taking a difference 



of 0.015 dex in abundance for a range of ~ 1.5 in the logio(W//l) 
of Fe lines. Then, wic\ ~ 0.25 and thus W2 ~ 2500. For the third 
condition, c\ = [FeI/H]-[FeH/H] and we adopt a maximum dif- 
ference of 0.015 dex. Then, W3C 2 ~ 0.25 and thus W3 ~ 1100 
dex -2 . Similarly, for W4 result W4 ~ 1100 dex -2 . In this estima- 
tion the weights Wi,...,W 4 resulted 18000 eV 2 /dex 2 , 2500, 1100 
dex -2 and 1 100 dex -2 , respectively, for a solution in which the 4 
conditions contribute equally with 0.25 to the x 1 function in the 
limit case of x 1 =1- In this example those solutions with^ 2 >1 
do not verify the four conditions. 

It is probably that the user have their own criteria adopt- 
ing the values of Wj,...,W4, instead of the example explained 
the previous paragraph. The user is free to modify the values 
of Wi,...,W4 (file fundpar.par) and this could result in more (or 
less) restrictive conditions. The code use this values to define 
X 1 and then search the minimum of the function. Then, the user 
should read the values of the slopes and metallicities in the out- 
put files to verify if the 4 conditions are satisfied. The values of 
Wj,...,W4 previously showed seems to verify in practice the re- 
queriments of minimization and verification of the 4 conditions. 
In the Table[TJwe show a sample of the file fundpar.par where the 
weights could be modified. Other parameters will be explained 
in the next sections. 

Following the definition, x 2 could be considered as 
a function that depends of the fundamental parameters 
X 2 =x 2 (T e ffJogg,[Fe/H],£;). If x 2 is not minimum, the algorithm 
should determine the next set of 4 possible values. These new 
variables are used in another iteration step (following the steps a 
to e) to derive a new model atmosphere, metallicity and finally 
a new value of x 2 ■ The algorithm that determine the next group 
of 4 parameters is the downhill method, explained in the next 
section. 

In the Table [2] we show a list of the input and output files 
used by FUNDPAR. The format of the input/output files is de- 
tailed in the file install.txt. There are two main directories (datain 
and dataout) containing the input and output files of the stars. 
The equivalent widths should be stored in separate files (one file 
by star), and the names of these files should be listed within an- 
other file called filenames.txt. The files atlas.par, batch.par and 
fundpar.par determine the value of some parameters used in the 
model calculation and abundance determination and will be ex- 
plained in the next sections. After the execution of FUNDPAR, 
there are three output files by star: the ATLAS model atmosphere 
of the solution and two output files directly from the MOOG 
abundance determination. The file output 1. screen contain infor- 
mation similar to the screen and output2. results list the final pa- 
rameters and their uncertainties. 



3. The downhill simplex method: minimization of \- 

In this section we briefly review the minimization procedure of 
X 2 as a function of 4 independ ent variables , using a Numerical 
Recipe's routine calle d amoeba (|Press 1 1992b. The downhill sim- 
plex method is due to lNelder & Mead d 19651) and requires only 
function evaluations, not derivatives. A simplex could be consid- 
ered as a geometrical figure of N+l vertices in a N-dimensional 
space (in our case, N=4). Taking any vertice as the origin, then 
the 4 other points define possible vector directions in the 4- 
dimensional volume. 

The downhill simplex method start with a group of N+l i.e. 
5 vertices rather than a single point or vertice. These vertices are 
desplaced in a characteristic length scale of the problem. In our 



case, x is initially calculated adopting displacements of 200 K, 
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tl 


200.00 


# Characteristic length scale of Teff (K) 


t2 


0.10 


# Characteristic length scale of Logg (dex) 


t3 


0.10 


# Characteristic length scale of [Fe/H] (dex) 


t4 


0.10 


# Characteristic length scale of Xita (km/s) 


wl = 


25000.00 


# Chi2 Weight factor 


w2 = 


2500.00 


# Chi2 Weight factor 


w3 = 


1111.00 


# Chi2 Weight factor 


w4 


1111.00 


# Chi2 Weight factor 



Table 2. List of the input and output files used by FUNDPAR. 



Filename 


Directory 


Comment 


Input Files: 


filenames.txt 
hd001.txt 


datain 
datain 


List of files containing the equivalent widths. 
Equivalent widths of the star HD 001 


hd999.txt 


datain 


Equivalent widths of the star HD 999 


Parameter Files: 


atlas. par 

batch.par 

fundpar.par 


atlas.launcher 

FUNDPAR 

FUNDPAR 


Parameters used by ATLAS 
Parameters used by MOOG 
Parameters used by FUNDPAR 


Output Files: 


hdOOl.mod 

hd001.outl.txt 

hd001.out2.txt 


dataout 
dataout 
dataout 


Model atmosphere of HD 001 from ATLAS 9 
Output file 1 from MOOG 
Output file 2 from MOOG 


output 1. screen 
output2.results 


dataout 
dataout 


Similar data to the screen 
Table of results for all the stars 



0. 1 dex, 0.1 dex and 0.1 km/s for the variables T e ff, logg, [Fe/H] 
and respectively. The characteristic lengths could be modi- 
fied in the file fundpar.par (see Table [TJ. FUNDPAR execute an 
initial calculation of x 2 at the vertices of the simplex, previous 
to the iteration process. Then, the downhill method takes a se- 
ries of steps, ussually moving the highest point of the simplex 

1. e. where x 2 is maximum. Succesive steps could be visualized 
as reflections, expansions and contractions of the 4-dimensional 
object. When the simplex found a valley, it contracts itself down 
the valley. An appropriate sequence of these steps will converge 
to the minimum of y 2 . 

As explained bv lPressI d 1992b . the termination criteria is us- 
sually delicate in the minimization process. The program require 
that the decrease in the function value in the terminating step 
be fractionally smaller that some tolerance (variable ftol within 
the amoeb a subroutine). The method explained exactly follows 
iNelde r & Mead ( 1965). However, in practice we found that for 
some stars the routine converge in a few iteration steps to a solu- 
tion with^K 2 »1. In such cases the code decrease the tolerance 
and continue the iteration process from the last point. We veri- 
fied that it is not necessary in this case to restart again all the pro- 
cess, because the effect of modify the tolerance determine only 
in which iteration the program stops. We also note that decrease 
the tolerance do not guarantee the convergence of the solution: 
we are only modifying the termination step or the termination 
criteria. 

Using the adopted values of Wi,...,W4, the program ussually 
takes less than ~200 iterations to reach a solution. The number 
of iterations is a known problem of the downhill method. If the 
code reach a solution with^ 2 >1, then the program restart the 
iteration process using the last so lution as a ne w initial condi- 
tion. The restart is recomended by Press] (1 19921) in the downhill 
method to eliminate a probable local minimum of the function. 
The code restart the iteration process only once. Also the user 



could manually restart again all the process using the last solu- 
tion, for instance, as new initial condition and FUNDPAR will 
try to search a solution with a smaller^ 2 . 



4. Derivation of the model atmospheres 

Together with the FUNDPAR code, we provide another indepen- 
dent program called atlas.launcher, which prepare and execute 
the Kurucz's ATLAS9 in LTE. The input of atlas.launcher are 
the fundamental parameters T e ff, logg, [Fe/H] and The out- 
put is an ATLAS solar-scaled model atmosphere corresponding 
to these parameters. The program use a parameter file called at- 
las.par that determine different parameters used by ATLAS, such 
as the mixing-length parameter (ussually taken as 1.25), and 
the ov ershooting weight parameter W, defined in ICastelli et al. I 
(119971) . A list of files used by ATLAS (grid of precalculated 
models, Rosseland and ODFs, and a sample of atlas. par) is pre- 
sented in the file install.txt. 

Then, the model atmosphere and the equivalent widths are 
introduced in the MOOG program. We instructed thi s code with 
the solar abunda nces of Grev esse & Sauvall d!998l) instead of 
the original from I Anders & Grevesse I d 1989h . except for Fe for 
which we adopted [Fe/H] - 7.47 dex. The N EWODF opacities 
use these abundance values (ICastelli & Kurucz 1 12003) . MOOG 
use a parameter file (batch.par) where some optional could be 
modified, such as the molecular equilibrium and the van der 
Waals line damping options. 

Literature authors have different preferences in the choice of 
ATLAS and MOOG parameters. The user is free to modify the 
parameters under their own criteria, and then FUNDPAR will 



4 For a complete list of options, see the manual of the MOOG pro- 
gram, http ://verdi . as . utexas . edu/moog. html 
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Table 3. Parameters adopted in different executions A, 
sample stars. (See text for details) 



,E for 5 



Run 


ML 


W 


D 


A 


1.25 


1 


2 


B 


1.25 





2 


C 


1.32 


1 


2 


D 


1.25 


1 


1 


E 


1.25 


1 






found the corresponding solution. To give an idea of the sensiv- 
ity of the results, we rederived the fundamental parameters of 5 
sample stars (HD 106252, HD 177830, HD 190228, HD 195019 
and HD 202206) using different combination of the parameters. 
In the Table [3] we show the values of the parameters adopted 
in different calculations A,...,E. The parameters showed are the 
mixing-length parameter (ML), the overshooting parameter (W) 
and the line damping option (D). We start the execution A adopt- 
ing ML=1.25, W=l and D=2. In the execution B we switched off 
the overshooting, while in the execution C we adopted a slightly 
higher ML. In the executions D and E we used different options 
for the line damping. 

The variation in T e ff> log g, [Fe/H] and £ (adopting different 
ML, W and D) is different from star to star. For instance, vary- 
ing W some stars increase the T e ff while others decrease the 
T e ff. Probably the variation depends of the fundamental param- 
eters of the stars. Then, in the Figure Q] we show the difference 
(called A) between the parameters derived in the A execution 
and the parameters derived in the executions B, C, D and E. The 
differences are plotted vs. T e ff only f° r the 5 stars previously 
mentioned. All differences are derived with respect to the A ex- 
ecution. The symbols used in the panels are filled circles (results 
B-A), diagonal crosses (results C-A), empty circles (results D- 
A) and squares (results E-A). The panels seems to show a slight 
tendence with T e ff, Switching on/off the overshooting (B execu- 
tion), AT e yy, Alog g and A[Fe/H] seems to show a tendence with 
T e ff- Increasing ML from 1.25 to 1.32 (C execution), we find no 
clear tendences and the differences are small. Using the damp- 
ing option 1 instead of 2 (D execution), T e ff, logg and [Fe/H] 
increased an average of ~48 K, ~0.14 dex and ~0.05 dex, re- 
spectively. Using the damping option (E execution), T e ff> l°gg> 
[Fe/H] and £ increased an average of ~70 K, ~0.20 dex, ~0.06 
dex and ~0. 1 km/s, respectively. We caution that these prelim- 
inar differences have been derived using only 5 stars. A higher 
number of stars is needed to properly determine how the param- 
eters ML, W and D modify the derived fundamental parameters. 



5. Uncertainty estimation of the derived parameters 

The unce rtainties in the atmospheric p arameters were estimated 
following lGonzalez & Vanture I (fl99 8jl The uncertainty in £ was 
determined from the standard deviation in the slope of the least- 
squares fit of abundance vs. reduced equivalent width. Then, the 
dispersion in effective temperature was determined from the un- 
certainty in the slope of the least-squares fit of abundance vs. 
excitation potential, in addition to the uncertainty in the slope 
due to the uncertainty in The uncertainty in the Fe abundance 
was derived combining the uncertainties in T e ff, £, and the scat- 
ter of the individual Fel abundances (standard deviation of the 
mean), all added in quadrature. In calculating the uncertainty in 
logg, we include the contribution from the uncertainty in T e ff in 
addition to the scatter in the Fell line abundances. 



We derive another estimation of the uncertainty of the solu- 
tion using x 2 - As we explained previously, we selected the values 
of the weights such that solutions with x 1 > 1 do not verify the 
4 conditions. Then, we adopt the size of the region^ 2 =1 as an- 
other estimation of the uncertainty of the solution. FUNDPAR 
store a record of the points with^ 2 <1 and use them to estimate 
the size of the region. The range of the values of T e ff, logg, 
[Fe/H] and £ are showed in the results with the number n of so- 
lutions with^ 2 < 1 . These dispersion values should be taken with 
caution if n is small. This kind of uncertainty is com parable to 
those derived by the criteria of iGonzalez & Vanture I d!998l) . In 
the next section we present a histogram comparing both uncer- 
tainties for a group of exoplanet host stars. The user is free to 
select between them, or the maximum value, for instance. 

FUNDPAR include the uncertainty estimation of the param- 
eters and they are showed in the ouput files. The amount of the 
uncertainty depend on many factors, such as the number of lines 
involved (which is ussually low for Fell), the measured equiv- 
alent widths and the laboratory data of the spectral lines. The 
log gf and excitation potential of the lines are important because 
many derived values (Fe abundances, the slopes of abundance 
vs. equivalent widths and vs. excitation potential, etc.) depends 
on these quantities. We note that the inclusion of a line which 
(for one reason or another) result in an abundance very different 
from the average, could modify the results and/or increase the 
uncertainty of the parameters. These lines should be eliminated 
from the list of measured equivalent widths. The abundance of 
individual lines are showed in the file out2.txt (work.dir direc- 
tory). 



6. Comparison with literature 

There are small differences in the paramete r s der ived using this 
method by different authors. ISantos et al. I d2000h derived sur- 
face gravities systematically highe r than the ones obtained by 
other authors (e.g. Gonzalez et al. 1 12001 by 0.15 dex. To solve 
this problem, ISantos et al. adopted new log gf values 

for the iron lines. Then the authors found a small average dif- 
ference of +25 K in tem perature compared to the studies of 
IGonzalez et al. I d200lh and lLaws et al. I (12003h (57 stars in com- 
mon). Also their log gs are +0.05 dex (on average) above and the 
avera ge differences in metallicity are between -0.10 and +0.10 
dex. lAmmler-von Eiff et al.1 (120091) derived the fundamental pa- 
rameters of host stars with transiting planets, using the 2002 
version of the code MOOG and ATLAS model atmospheres. 
However f or the TrES and HAT objec ts, the abundan c e det ermi- 
nations of lSozzetti etal.l d2007l) and ISozzetti et al. I d2009l) are 
systematically lower than those derived in their work by ~0.10 
dex, using essentially the same method. The authors also men- 
tion differences in temperature (~100 K) and gravity (~0.15 dex) 
due to possible sistematic tendences between the works, and they 
note that the origin of the discrepancies in abundances is com- 
monly unidentified. 

We compare the p arameters derived using FUNDPAR and 
those from literature dGonzalezlfl997l 1998; Gonzalez et al. 



1 19991: [Santos et al. I ~~l2000t IGonzalez & Laws I 1200 
Gonz alez et al. I 1200 U lLaws & Gonzalez! 120011: lLaws et al. 
2003). We selected these works in particular because they 
present the Fe equivalent widths: we use the same values with 
FUNDPAR. In the Table|4]we show the T eff , logg, [Fel/H] and 
£ obtained. Some stars in the table are listed twice (identified 
with a subscript) because we have different equivalent width 
sources from literature. The sources of the equivalent widths are 
listed in the Table g] as R1,...,R8. 
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Fig. 1. Difference between the parameters derived in the A execution and the parameters derived in the executions B, C, D and 
E. The differences are plotted vs. T e ff on ly f° r 5 stars. For instance, in the T e /f panel we ploted AT e yy = T e yy(B)-T e //(A) with 
filled circles. The symbols used in the panels are filled circles (differences B-A), diagonal crosses (differences C-A), empty circles 
(differences D-A) and squares (differences E-A). 



We show in the Figure [2] an example of abundance vs. exci- 
tation potential and abundance vs. logw(W/A) (logarithm of the 
reduced equivalent width) for HD 106252. In the Figure [3] we 
present the difference between FUNDPAR parameters and litera- 
ture (y axis, A=FUNDPAR-literature) vs. literature (x axis). The 
empty point in the Figure [3] correspond to the star HD 38529, 
for which their slopes, metallicities and x 1 values indicate that 
this object should be taken with caution. The derived parame- 
ters of this star are based only on 3 Fell lines (and 24 Fel lines). 
In literature, some values have been rounded in T e ff> log g an d 
£ to within 50 K or 10 K, 0.05 dex and 0.1 km/s, respectivel y 
(see, for example. fGonzalez lll997lll998tlGonzalez et al. Ill999l) . 
There is a good agreement between FUNDPAR parameters and 
previous works from literature within the errors, which is logic 
taking into account that FUNDPAR use a very similar method. 
However the Figure shows that there is a dispersion in the val- 
ues of the parameters (particularly logg and £) and probably a 
sistematic tendence for the metallicity (~0.01 dex below litera- 
ture values, see next discussion). The values of AT e // and Alog g 
also seems to slightly decrease with T e /f and log g, respectively. 
The median of the differences for the fundamental parameters 
compared to literature are 24 K, 0.06 dex, 0.03 dex and 0.08 
km/s, corresponding to T e ff> l°gg, [Fe/H] and respectively. 
The higher differences in the parameters are 118 K (16 Cyg B), 
0.30 dex (HD 27442), 0.16 dex (16 Cyg B) and 0.31 km/s (47 
Uma), corresponding to T e ff, logg, [Fe/H] and respectively. 



Now we discuss the possible origin of the dispersions 
and probable tendences observed in the Figure [3] We tested 
FUNDPAR modifying significatively the values of the weights 
Wj,...,W4 in the function x 2 an d verified that the difference 
in [Fe/H] with literature and the dispersion in the other pa- 
rameters changes very slightly. Then the weights Wi,...,W4 do 
not seem to be the cause. We use a method similar to liter- 
ature, however it is not totally identi cal. FUNDPAR use dif- 
ferent Kurucz model atmospheres (Caste lli & Kurucz 1 12003h 
than those used in literature (most of them are previous to the 
creation of the ODFNEW models). The code use model at- 
mospheres derived by ATLAS using ODFNEW opacities and 
solar abundances from iGrevesse & Sauval I d!998l) instead of 
lAnders & Grevessel d 1989b . The ne w models presen t differences 
compared to older Kurucz models dKurucz II 19901) such as the 
solar abundances, the replacement of TiO and H2O molecular 
lines, some HI quasi-molecular absorptions, etc. taken into ac- 
count in the NEWODF opacities. Preliminar improvements are 
the U-B and u-b color indices for T e ff <6750 K, all color in- 
dices for cooler stars, and the better modeling for the u pper lay- 
ers of cool and giant stars (ICastelli & Kuruczll2003l) . In this 
example, the model atmospheres are computed with convec- 
tion (mixing-lenght parameter= 1 .25), overshooting (W=l) and 
damping in the Unsold approximation b ut multiplied by a factor 
as suggested by Bl ackwell et al. I d 19951) . The use of the fourth 
condition within the function^ 2 is (possibly) another difference 
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Fig. 2. Plots of abundance vs. excitation potential and abundance vs. logw(W/A) (logarithm of the reduced equivalent width) for HD 
106252. 



with previous studies. Literature works surely take this into ac- 
count, however it is not totally clear for us how. Finally, we use 
the MOOG 2009 versiorQ of the program instead of the 2002 
version used in literature, although we expect almost the same 
abundance values from both versions. 

These differences, at least in part, produce the slightly 
disimilar values showed in the Figure [3] which suggest that 
FUNDPAR use probably a different metallicity scale than used 
in literature. A complete comparison require the exact knowl- 
edge of all involved details used in the literature calculation. 
Our intention is to clearly present all the asumptions used in 
FUNDPAR, in the model atmospheres and within the code. 

In The Figure [4] we show the histogram distributions of the 
uncertainties derived in T e ff, logg, [Fe/H] and ^. The densely 
and slightly sh aded histograms corre s pond to uncertaities de- 
rived following Gon zalez & Vanture ] (fl998h and using the x 
function, respectively. Some distributions present a peak in a 
common uncertainty value, such as ~0.05 dex in the [Fe/H] dis- 
tribution and ~30 K in the distribution of effective temperature. 
We see that the errors derived using both methods are compara- 
ble. 

The metallicities presented in the Table [4] correspond to a 
group of exoplanet host stars. The median of the group is 0.17 
dex with a dispersion of 0.22 dex. In the Figure [5] we present 
the histogram of the metallicity distribution. Then, as an exam- 



0.3 




Fig. 5. Metallicity distribution of the sample stars derived by 
FUNDPAR. 



http ://verdi . as .utexas . edu/moog. html 



pie of practical use of FUNDPAR, we verified the metal-rich 
nature of main sequence stars with low mass c ompanions , a fact 
known from the li terature (see, for example, iGonzalezI Il997t 
ISantos et al. Il2000h . 
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Fig. 3. Difference between FUNDPAR parameters and literature (y axis) vs. literature (x axis). The panels shows separately the 
effective temperature, gravity, metallicity and Filled and empty points correspond to solutions with \ 2 < 1 and x 2 > 1 , respectively. 



7. Concluding remarks 

We have implemented a fortran algorithm available from the 
web that estimate fundamental parameters of solar type stars, 
requiring only the measure of Fe equivalent widths. The fi- 
nal solution should verify the three conditions of the standard 
method: [FeI/H]=[FeII/H] (i.e. ionization equilibrium), indepen- 
dence of the metallicity with the excitation potential (i.e. exci- 
tation equilibrium) and with respect to the equivalent widths. 
We also add another condition: the input metallicity used in the 
model atmosphere should be similar to the resulting metallicity 
from the equivalent widths. We taken into account these condi- 
tions in one variable called^ 2 , adopting an expression which in- 
clude the weights Wj,...,W4. FUNDPAR use Kurucz model atmo - 
spheres with the NEWODF op acities dCastelli & Kurucz 1120031) . 
solar-scaled abundances from iGrevesse & Sauvall (1998) and 
the 2009 version of the MOOG program. Different details could 
be selected, such as the mixing-lenght parameter, the overshoot- 
ing and the damping of the lines, for instance. We have planed a 
new version that include the option of use the WIDTH9 program 
instead of MOOG deriving abundances from equivalent widths. 

The code include the derivation of the uncertainty in the 4 
parameters following the criteria of Gonzalez & Vanture (1998) 
and another uncertainty estimation using the x 1 function. We 
verified the metal-rich nature of a group of exoplanet host stars. 
The parameters derived with FUNDPAR are in agreement with 
previous works in literature. 
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Fig. 4. Histogram distributions of the uncertainties in T e ff, logg, [Fe/H] and £ derived with FUNDPAR. The densely and slightly 
shaded histograms correspond to uncertaities derived following lGonzalez & Vanture l (fl998h and using the^ 2 function, respectively. 
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Table 4. Fu ndamental parame ters a nd dispersions derived with FU NDPAR. Referen c es for the e quivalent widths (la st col umn): Rl: Gonzalez et al. 
dl999h. R2:lGonzalez dl998h. R3:lLaws & Gonzalez I d200lh . R4: iGonzalez et ail d200lh . R5 : Laws et al. I d2003l) . R6: ISantos et al. I d2000h . R7: 
iGonzalezI (2000), R8: Gonzalez] d!997l) 



Star 


T eff 


logg 


[Fe/H] 




eTeff 


elogg 


e[Fe/H] 




Ref 




[K] 


[dex] 


[dex] 


[km/s] 


[K] 


[dex] 


[dex] 


[km/s] 




14 Her 


5294 


4.33 


0.50 


0.71 


40.60 


0.19 


0.07 


0.05 


Rl 


16 Cyg Ai 


5704 


4.29 


-0.01 


1.32 


90.75 


0.19 


0.11 


0.10 


R2 


16 Cyg A 2 


5749 


4.23 


0.07 


1.02 


15.58 


0.04 


0.06 


0.03 


R3 


16CygBi 


5632 


4.36 


-0.05 


1.21 


83.68 


0.25 


0.10 


0.09 


R2 


16 Cyg B 2 


5690 


4.30 


0.05 


0.91 


14.16 


0.04 


0.06 


0.03 


R3 


47 Uma 


5720 


4.26 


-0.12 


1.31 


86.63 


0.15 


0.10 


0.08 


R2 


51 Pegi 


5681 


4.37 


0.11 


1.21 


62.21 


0.11 


0.08 


0.08 


R2 


51 Peg 2 


5797 


4.39 


0.22 


1.13 


16.33 


0.04 


0.04 


0.03 


R4 


70 Vir 


5481 


3.98 


-0.08 


1.10 


92.16 


0.15 


0.09 


0.07 


R2 


BD-10 3166 


5275 


4.36 


0.32 


0.69 


42.70 


0.05 


0.09 


0.05 


R4 


HD 106252 


5844 


4.50 


-0.06 


0.98 


12.23 


0.12 


0.05 


0.03 


R5 


Tir\ 1 r\ r\*i 

HD 10697 


5600 


3.96 


0.10 


1.20 


28.41 


0.03 


0.05 


0.04 


R4 


HD 108147 


6279 


4.59 


0.19 


1.11 


41.51 


0.12 


0.11 


0.09 


R5 


HD H4783 


5100 


4.48 


0.08 


0.87 


15.70 


0.10 


0.07 


0.04 


R5 


HD H7176 


5487 


4.09 


-0.06 


1.00 


10.86 


0.10 


0.06 


0.02 


R5 


HD I237i 


5537 


4.50 


0.22 


1.19 


30.76 


0.15 


0.07 


0.05 


R4 


HD 12372 


5468 


4.50 


0.05 


1.38 


23.27 


0.18 


0.07 


0.04 


R6 


HD 12661 1 


5670 


4.41 


0.38 


0.91 


19.18 


0.06 


0.05 


0.03 


R4 


HD 126612 


5719 


4.48 


0.38 


0.84 


24.28 


0.04 


0.05 


0.04 


R4 


HD 13445 


5080 


4.43 


-0.30 


0.65 


19.69 


0.23 


0.07 


0.05 


R6 


HD 134987 


5719 


4.28 


0.33 


1.09 


22.44 


0.04 


0.07 


0.04 


R4 


HD 136118 


6153 


4.37 


-0.09 


2.00 


21.82 


0.15 


0.08 


0.11 


R5 


HD 141937 


5832 


4.50 


0.12 


0.97 


12.63 


0.10 


0.05 


0.03 


R5 


HD 160691 


5788 


4.49 


0.24 


1.12 


22.03 


0.08 


0.06 


0.04 


R5 


HD 16141 1 


5768 


4.24 


0.17 


1.10 


22.07 


0.04 


0.04 


0.04 


R4 


HD 16141 2 


5751 


4.28 


0.16 


1.00 


12.18 


0.05 


0.05 


0.03 


R5 


HD 168443 


5549 


4.13 


0.08 


1.04 


17.06 


0.02 


0.06 


0.03 


R4 


HD 168746 


5566 


4.43 


-0.09 


0.94 


10.51 


0.04 


0.05 


0.03 


R5 


HD 169830] 


6268 


4.10 


0.20 


1.29 


22.14 


0.18 


0.05 


0.03 


R6 


iir\ 1 r\ ^1 f\ 

HD 1698302 


6283 


4.23 


0.16 


1.27 


24.27 


0.08 


0.06 


0.06 


R5 


HD 177830 


4801 


3.37 


0.38 


0.95 


55.16 


0.07 


0.10 


0.04 


R4 


HD 187123 


5792 


4.39 


0.14 


1.04 


26.26 


0.08 


0.04 


0.04 


Rl 


HD 190228 


5259 


3.64 


-0.32 


1.13 


11.79 


0.11 


0.05 


0.02 


R5 


HD 192263 


4905 


4.38 


-0.02 


0.73 


19.76 


0.18 


0.08 


0.05 


R4 


HD 195019a 


5731 


4.17 


0.00 


1.20 


13.09 


0.13 


0.05 


0.03 


R5 


HD 19994.2 


6147 


4.35 


0.12 


2.03 


23.94 


0.10 


0.08 


0.09 


R5 


HD 202206 1 


5635 


4.55 


0.28 


0.88 


24.27 


0.14 


0.06 


0.04 


R6 


HD 2022062 


5713 


4.45 


0.30 


1.07 


14.42 


0.04 


0.06 


0.03 


R5 


HD 209458 


6065 


4.43 


0.02 


1.14 


18.56 


0.13 


0.06 


0.04 


R4 


HD 210277 


5517 


4.34 


0.24 


0.84 


26.56 


0.12 


0.05 


0.04 


Rl 


HD 217107 


5602 


4.42 


0.42 


0.87 


17.64 


0.03 


0.05 


0.03 


R4 


HD 22049 


5083 


4.47 


-0.14 


0.84 


14.29 


0.15 


0.06 


0.04 


R5 


HD 222582 


5719 


4.25 


0.00 


0.99 


18.90 


0.03 


0.04 


0.04 


R4 


HD 27442 


4859 


3.60 


0.38 


1.21 


45.85 


0.22 


0.15 


0.06 


R5 


HD 28185 


5637 


4.54 


0.21 


0.96 


13.79 


0.04 


0.06 


0.03 


R5 


HD 33636 


5877 


4.27 


-0.13 


0.95 


14.81 


0.10 


0.06 


0.04 


R5 


HD 37124i 


5579 


4.60 


-0.39 


0.74 


19.83 


0.11 


0.09 


0.08 


R5 


HD 37124 2 


5495 


4.46 


-0.42 


0.71 


22.54 


0.04 


0.05 


0.07 


R4 


HD 38529 


5600 


3.82 


0.34 


1.25 


38.21 


0.05 


0.06 


0.03 


R4 


n U 4ZUJ) 


JOO 1 




fi A*X 


1 MG 
1 .Uo 


z.y.yj 


M MA 
U.UO 


M 1 ^ 


M MA 
U.UO 


Kj 


HD 4208 


5592 


4.33 


-0.29 


0.87 


14.74 


0.15 


0.06 


0.04 


R5 


HD 46375 1 


5210 


4.37 


0.24 


0.65 


35.26 


0.10 


0.06 


0.04 


R4 


HD 46375 2 


5252 


4.50 


0.26 


0.68 


21.15 


0.09 


0.10 


0.06 


R5 


HD 50554 


5978 


4.51 


-0.01 


1.10 


11.93 


0.10 


0.06 


0.03 


R5 


HD 52265 1 


6102 


4.28 


0.24 


1.20 


17.71 


0.06 


0.05 


0.03 


R4 


HD 52265 2 


6163 


4.37 


0.27 


1.21 


19.51 


0.10 


0.03 


0.03 


R4 


HD 52265 3 


5993 


4.21 


0.17 


1.24 


25.17 


0.18 


0.06 


0.04 


R6 


HD 6434 


5715 


4.32 


-0.54 


0.70 


23.11 


0.16 


0.11 


0.10 


R5 
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Table 4. Continued. 



Star 


J eff 


logg 


[Fe/H] 




eTeff 


elog g 


e[Fe/H] 




Ref 




[K] 


[dex] 


[dex] 


[km/s] 


[K] 


[dex] 


[dex] 


[km/s] 




HD 68988 


5936 


4.49 


0.32 


1.19 


19.40 


0.02 


0.09 


0.05 


R5 


HD 75289 1 


6138 


4.52 


0.27 


1.50 


20.96 


0.10 


0.07 


0.05 


R7 


HD 75289 2 


6070 


4.44 


0.23 


1.42 


30.93 


0.16 


0.06 


0.05 


R6 


HD 75332 


6277 


4.47 


0.22 


1.03 


24.73 


0.11 


0.05 


0.04 


R4 


HD 82943 1 


5950 


4.45 


0.28 


1.10 


25.37 


0.07 


0.05 


0.03 


R6 


HD 82943 2 


5987 


4.46 


0.24 


1.11 


17.85 


0.04 


0.05 


0.03 


R5 


HD 83443 1 


5389 


4.31 


0.33 


1.00 


27.40 


0.18 


0.08 


0.04 


R6 


HD 8 3443 2 


5427 


4.40 


0.34 


0.99 


22.07 


0.13 


0.07 


0.05 


R5 


HD 8574 


6029 


4.34 


0.00 


1.19 


13.29 


0.16 


0.05 


0.03 


R5 


HD 89744 


6315 


4.16 


0.27 


1.67 


28.06 


0.07 


0.04 


0.05 


R4 


HD 92788 


5755 


4.43 


0.29 


1.02 


26.37 


0.04 


0.05 


0.04 


R4 


HD 95128 


5832 


4.31 


0.02 


1.08 


13.83 


0.10 


0.05 


0.02 


R5 


HR 810 


6072 


4.48 


0.15 


1.23 


27.89 


0.09 


0.07 


0.08 


R4 


v Andi 


6165 


4.21 


0.13 


1.40 


50.44 


0.03 


0.06 


0.07 


R8 


v And2 


6150 


4.14 


0.09 


1.48 


50.48 


0.03 


0.07 


0.07 


R7 


p Cnci 


5201 


4.27 


0.32 


0.78 


99.87 


0.20 


0.11 


0.10 


R2 


p Cnc 2 


5206 


4.40 


0.41 


0.69 


27.42 


0.18 


0.07 


0.04 


R7 


pCrb 


5751 


4.15 


-0.24 


1.33 


28.85 


0.05 


0.04 


0.05 


R2 


T BOO] 


6517 


4.37 


0.30 


1.57 


60.99 


0.04 


0.06 


0.08 


R8 


T B002 


6415 


4.23 


0.30 


1.34 


49.61 


0.06 


0.08 


0.07 


R7 
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